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possible  to  construct  strong,  lightweight,  high  horsepower 
DC  motors.  This  has  occasioned  a  reassessment  of  electrome¬ 
chanical  actuators  as  alternatives  to  comparable  pneumatic 
and  hydraulic  systems  for  use  as  flight  control  actuators 
for  tactical  missiles.  A  dynamic  equivalent  circuit  model 
for  the  analysis  of  a  small  four  pole  brushless  DC  motor  fed 
by  a  transistorized  power  conditioner  utilizing  high  speed 
switching  power  transistors  as  final  elements  is  presented. 
The  influence  of  electronic  commutation  on  instantaneous 
dynamic  motor  performance  is  particularly  demonstrated  and 
good  correlation  between  computer  simulation  and  typical 
experimentally  obtained  performance  data  is  achieved.  The 
model  is  implemented  in  CSHP  language  and  features 
improvements  in  transistor  and  diode  models  as  well  as  a 
more  accurate  air  gap  flux  representation  over  previous 
work.  Hall  effect  sensor  rotor  position  feedback  is  simu¬ 
lated.  Both  constant  and  variable  air  gap  flux  is  modeled, 
and  the  variable  flux  model  treats  the  flux  as  a  fundamental 
and  one  harmonic. 
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I.  IHTBODCCTIQI 


Electronically  commutated  brushless  DC  aotors  are 
increasingly  becoming  practical  for  numerous  applications. 
The  advantages  of  these  aotors  over  conventional  brush-type 
£C  aotors  are  smaller  size  and  lighter  weight  per  horse¬ 
power.  The  increase  in  power  over  conventional  DC  motors  is 
realized  by  the  use  of  three  phase  circuitry  commonly  found 
in  synchronous  1C  aotors  which  results  in  additional  motor 
cur rent -magnetic  fluz  interaction  not  available  with  conven¬ 
tional  DC  motor  designs.  The  elimination  of  the  brush 
commutation  and  rotating  armature  allows  operation  at  higher 
speed  with  higher  currents  and  improves  the  thermal  charac¬ 
teristics  of  the  motcr  as  heat  generating  windings  will  be 
contained  in  the  stator  which  can  be  more  efficiently 
cooled. 

In  this  thesis,  an  equivalent  circuit  motor- power  condi¬ 
tioner  modeling  approach  is  used  to  predict  the  dynamic 
performance  of  a  typical  brushless  DC  motor  for  a  fin 
control  actuator  assembly.  The  interaction  of  the  power 
conditioner  final  drive  elements  with  the  stator  currents  is 
demonstrated  in  addition  to  the  rotor  magnetic  flux  -  stator 
current  interaction  and  interdependence.  Numerical  results 
of  steady  state  operation  are  compared  with  experimentally 
obtained  instantaneous  voltages  and  currents  in  amplitude 
and  in  profile.  After  demonstration  of  steady  state  opera¬ 
tion  cf  the  model,  start-up  and  reversal  transients  are 
predicted. 

The  model  developed  in  this  thesis  is  part  of  a  contin¬ 
uing  program  to  develop  a  comprehensive  simulation  for  a  fin 
control  actuating  system  which  has  been  the  subject  of 
previous  work  [Bef.  1  and  2]  and  is  subject  to  further 
development. 
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The  object  of  this  thesis  is  to  extend  previous  work  on 
brushless  DC  motors  [Bef .  1 ]  to  include  switching  character¬ 
istics  of  power  transistors  and  diodes  in  such  a  manner  as 
to  provide  a  coherent  model  of  power  conditioner  final  drive 
elements  and  a  brushless  DC  motor.  Additionally,  the  simu¬ 
lation  provides  steady  state  as  well  as  instantaneous 
dynamic  performance  during  start-up  and  sudden  reversal  of 
the  motor.  The  IBM  Continuous  System  Modeling  Program 
(CSMP)  was  chosen  primarily  because  of  the  inherent  attri¬ 
butes  of  CSMP  for  modeling  dynamic  physical  systems. 
Inasmuch  as  previous  model  development  utilized  this 
language,  some  continuity  of  development  is  maintained. 
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II.  SYSTEM  DESCB1PTI0H 


A.  SYSTEM  BLOCK  DIAGBAH 

The  simulation  diagram  for  the  model  is  depicted  in 
figure  2.1  and  demonstrates  the  interrelation  between  the 
magnetic  flux  model,  the  torque  coefficient,  and  the  gener¬ 
ated  lack  eif  voltage.  The  nonlinear  nature  of  the  gener¬ 
ated  voltage  and  torque  is  indicated  by  the  multipliers  in 
the  diagram  for  angular  speed  times  the  magnetic  flux, 
¥  (t)  0  (t)  ,  used  in  the  back  emf  calculation  and  the  motor 
current  times  the  magnetic  flux,  I(t)((t),  used  in  the 
generated  torque  calculation.  The  three  phase  electrical 
simulation  calculates  three  phase  currents,  which  are  multi¬ 
plied  by  the  torque  coefficient  and  normalized  flux  model 
for  each  phase  resulting  in  a  torque  due  to  the  current  in 
each  phase.  The  motcr  torque  is  the  sum  of  the  three  phase 
torques,  and  the  net  torque  is  the  motor  torque  less  the 
load  torgue.  Motor  rctation  is  then  calculated  by  applying 
the  net  torque  to  the  rotor  inertia  and  viscous  friction. 
The  integral  of  the  rotor  angular  speed  is  calculated  for 
angular  position,  which  is  fed  to  the  flux  and  transistor 
logic  models  to  control  the  generation  of  conduction 
currents  for  each  phase.  The  transistor  model  and  magnetic 
flux  model  are  discussed  in  Chapter  Three  and  Chapter  Four, 
respectively. 


POWER  CONDITIONER 
LOGIC  CONTROLLER 


Figure  2.1  Motor  and  Power  Conditioner 
Simulation  Diagram 


The  turn-off  time  for  a  simple  switch  is  characterized 
by  the  following  equations:  [Ref.  4] 
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just  prior  to 

turn-off,  and 

I8,off  is  the  base  drive  current  immediately  following  turn¬ 
off  initiation.  An  inspection  of  equations  3.7  and  3.8 
reveals  that  the  turn-off  time  is  affected  by  twc  mecha¬ 
nisms,  one  of  base  drive  current  and  one  of  collector  or 
emitter  current  initial  condition.  Reducing  the  initial 
condition  of  collector  or  emitter  current  or  increasing  the 
base  drive  current  results  in  a  proportional  decrease  in 
turn-off  time. 


One 

of 
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difficulties  of  overdrive 

is  storage  time 

which  is  a 

function  of  the  excess  base 

current. 
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alt  +  oJot 

(3.  11) 
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and  0(u  corresponds  to  with  collector  and  emitter  termi¬ 
nals  reversed,  and  a)*  is  the  -3JB  alpha  cutoff  frequency 
with  collector  and  emitter  terminals  reversed.  [Ref.  4  and 

5]- 

It  is  most  desirable  to  minimize  the  storage  time  anl 
this  can  be  accomplished  by  overdriving  during  turn-cn, 
reducing  the  base  drive  to  hold  the  transistor  near 
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charge  is  swept  out  of  the  collector  and  emitter.  The 
minority  carrier  density  at  the  collector  is  again  small, 
and  the  charge  gradient  is  decreasing  rapidly.  As  soon  as 
minority  carrier  density  reaches  zero  at  the  emitter,  it  may 
also  become  reverse  biased.  This  turn-off  process  is  one  of 
the  important  limitations  on  the  switching  response  of 
transistors. 

The  turn-on  time  for  a  simple  switch  is  approximated  by 
the  following  equations:  [Ref.  4] 
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and  CJ^  is  the  -3dB  alpha  cutoff  freguency. 

The  effects  of  overdrive  and  reverse  drive  on  transistor 
switching  performance  are  demonstrated  in  Figure  3.2.  If 
not  overdriven,  it  takes  approximately  three  time  constants 
(Tcc  or  Tee)  to  reach  saturation.  A  significant  improvement 
in  rise  time  is  realizeable  by  the  use  of  base  overdrive. 
If  a  transistor  is  overdriven,  its  effective  gain  is 
decreased,  but  this  has  the  effect  of  slightly  reducing  the 
time  constant  while,  more  significantly,  attempting  to  drive 
the  transistor  to  higher  steady  state  current — resulting  in 
a  greater  initial  slope.  The  combination  of  the  two  effects 
results  in  rise  times  on  the  order  of  .7  time  constants  to 
reach  saturation. 


tase  current.  Eventually  the  charge  distribution  will 
bee o ae  Bore  unifora  with  a  constant  gradient  of  charge 
between  the  emitter  and  collector.  This  is  the  noraal 
active  transistor  region.  If  the  load  line  of  the  circuit 
is  such  that  the  transistor  is  allowed  to  pass  through  the 
active  region  and  enter  the  saturation  region,  the  collector 
current  reaches  a  maximum  value  and  cannot  increase.  The 
collector  to  base  voltage  is  then  very  snail,  but  the  base 
current  can  be  increased  further  so  that  the  minority 
carrier  distribution  in  the  base  region  has  essentially  zero 
gradient  fron  the  eaitter  to  collector.  This  creates  a 
large  excess  of  linority  carriers  at  the  collector  junction 
which  in  turn  causes  the  collector  to  be  forward  biased. 
This  is  a  state  in  which  there  is  more  charge  at  the 
collector  than  the  collector  can  collect.  In  this  condi¬ 
tion,  the  collector  to  emitter  voltage  is  smaller  than 
either  of  the  juncticn  voltages. 

The  amount  of  charge  stored  in  the  base  region  in  this 
condition  is  very  large.  additionally,  since  the  collector 
junction  is  forward  biased,  excess  minority  carriers  will 
cross  the  junction  from  the  base  region  and,  if  the  minority 
carrier  lifetime  is  sufficiently  long,  be  stored  in  the 
collector  region  as  well  as  in  the  base  region.  This  stored 
minority-carrier  charge  has  a  deleterious  effect  on  the 
turn-eff  time  of  the  transistor.  In  order  to  turn  off  the 
transistor,  the  minority  carrier  density  at  the  collector 
junction  must  be  reduced  to  zero  for  the  collector  voltage 
to  start  to  recover.  When  turn-off  begins,  the  charge 
gradient  is  smaller  at  the  emitter  than  at  the  collector, 
resulting  in  current  flow  into  the  base  lead  and  minority 
carriers  disappearing  from  the  base  region. 

The  collector  current  remains  unchanged  until  the 
storage  period  is  over  and  collector  voltage  remains  small. 
As  the  transistor  enters  the  active  region,  the  remaining 


B.  SI1TCHIHG  TBiMSISTOB  DYI1HICS 

Behavior  of  transistors  in  the  cutoff  and  saturation 
regions  differs  somewhat  froa  the  aore  typically  understood 
snail  signal  aodels.  Soae  discussion  of  switching  charac¬ 
teristics  and  the  resulting  nodel  is  in  order.  This  discus¬ 
sion  assuaes  a  generic  transistor  and  is  egually  correct  for 
PNP  and  HP!)  devices,  as  the  discussion  describes  primarily 
base  region  effects.  It  might  be  necessary  to  keep  in  mind 
that  electrons  are  majority  carriers  for  n-type  base  regions 
and  holes  are  majority  carriers  for  p-type  regions.  This 
discussion  of  transistor  behavior  follows  that  found  in 
£Hef.  4]. 

When  a  transistor  is  operating  in  the  cutoff  region,  the 
emitter  and  collector  are  reverse  biased,  base  region 
minority  carrier  concentration  is  practically  zero  at  the 
junctions,  and  each  junction  is  drawing  some  fraction  of 
leakage  current  Ico.  The  majority  carriers  which  are 
collected  at  the  junctions  are  thermally  generated  in  the 
base  region,  and  a  (small)  minority  carrier  current  will 
flow  cut  of  the  base  terminal  to  keep  the  base  region 
neutral.  The  amount  of  minority  carrier  charge  stored  in 
the  base  region  is  negligibly  small. 

To  turn  the  transistor  on,  charge  carriers  must  be 
injected  at  the  emitter  junction  so  that  there  will  be  a 
large  charge  gradient  toward  the  collector.  This  means  that 
there  must  be  minority  carriers  stored  in  the  base  region 
and  there  must  be  a  corresponding  change  in  majority 
carriers  to  neutralize  the  stored  charge.  Thus  current 
flows  into  the  base  region  through  the  base  terminal  as  the 
storage  takes  place.  If  the  base  current  is  suddenly 
increased,  there  will  be  a  large  charge  distribution 
gradiert  near  the  emitter  and  a  small  gradient  near  the 
collector,  the  difference  in  slope  being  proportional  to  the 
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TRANSISTOR  SWITCHING  LOGIC 


III.  ELECTBOHIC  COHHDTATIOH 


A.  TBBEI  PHASE  BHIDGI  CIBCOITBI 

The  three-phase,  four  pole  brushless  DC  notor  is 
depicted  in  Figure  2.2.  The  transistor  and  diode  configura¬ 
tion  is  very  siailar  in  appearance  to  a  three  phase  full 
nave  bridge  rectifier  circuit  and  acts  in  a  very  similar 
fashion.  The  composite  flux  shown  in  Figure  4.6  is  a  demon¬ 
stration  of  this  similarity.  The  composite  flux  actually 
corresponds  to  the  peaks  of  the  back  emf  voltage  waveform  in 
Figure  4.5  as  if  it  were  rectified  in  a  three  phase  full 
wave  bridge  rectifier.  This  figure  was  generated  by  taking 
the  normalized  flux  model  for  the  two  conducting  phases 
during  each  commutation  interval,  adding  them  together, 
multiplying  by  the  torque  coefficient  and  converting  to 
webers  (the  torque  coefficient  has  the  same  basic  units  of 
force-length/amp) .  Actually,  the  diodes  may  be  used  as  a 
full  wave  bridge  rectifier  for  some  applications  (e.g.  to 
recharge  batteries)  when  the  motor  is  driven  by  the  lead 
torque  and  acts  as  a  generator  [Hef.  3]. 

Transistors  Q1  through  Q6  control  the  current  switching 
in  accordance  with  the  timing  diagram  in  Figure  3. 1,  where  a 
logical  "I"  corresponds  to  the  transistor  being  in  the 
conduction  mode.  This  diagram  was  derived  for  counterclock¬ 
wise  rotation  (positive  direction)  as  the  current  will  be 
positive  in  each  winding  as  the  north  pole  of  the  rotor  is 
positioned  in  its  field  (positive  flux)  and  the  current  will 
be  negative  when  the  south  pole  of  the  rotor  is  positioned 
in  its  field  (negative  flux).  A  similar  timing  diagram  can 
be  made  by  considering  the  logic  for  clockwise  rotation 
(negative  direction)  in  Table  II. 
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listings  are  contained  in  Appendix  B.  The  first  of  these 
nodels  continues  to  siaolate  the  transistors  without  the 
protecting  diodes  (which  are  shown  in  Figure  2.2)  in  order 
to  deaonstrate  the  inductive  voltage  transients  produced  by 
opening  one  of  the  stator  winding  circuits  with  a  transistor 
switch.  The  next  ncdel  adds  the  daaper  diodes  to  the  power 
conditioner  to  fully  simulate  the  power  conditioner  final 
drive  eleaents.  The  final  aodel  takes  all  the  coaplications 
of  the  preceeding  ones  and  depicts  the  transient  behavior  of 
the  systea  during  reversal  ("plugging")  of  the  rotating 
■otor.  The  results  of  conputer  runs  using  this  mcdel  are 
discussed  in  Chapter  Five.  No  attempt  has  been  aade  to 
reduce  the  order  of  the  aodel  in  these  simulations  since  the 
desire  was  to  observe  both  the  fast  and  slow  time  constant 
effects.  The  numerical  integration  of  the  differential  equa¬ 
tions  of  such  a  "stiff"  system  to  evaluate  the  steady  state 
behavior  can  be  somewhat  time  consuming.  Transistor  time 
constants  are  user  adjustable  as  are  saturation  and  cutoff 
parameters.  The  aotcr  parameters  used  for  the  model  devel¬ 
opment  are  typical  fer  the  brushless  motors  now  available 
for  this  type  of  application,  but  are  not  representative  of 
any  particular  motor.  As  such,  all  motor  parameters  are 
user  adjustable  to  conform  to  any  desired  motor  performance 
that  would  meet  design  objectives. 
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aotor  reversal  cannot  be  accoaplished  by  reversing  the 
polarity  of  the  pover  supply  as  in  ££e£.  1]  but  aust  instead 
be  accoaplished  as  in  a  typical  power  conditioner  by  proper 
reverse  sequencing  of  the  connotation,  shown  in  Table  II, 
and  discussed  in  Chapter  5. 

Several  nodels  were  produced  whereby  each  nodel  repre¬ 
sented  a  coapletion  cf  a  particular  aspect  of  the  systen 
such  as  aotor  air  gap  flux  and  power  conditioner  transistor 
switching  characteristics.  The  early  nodels  were  developed 
froa  siaple  ideas  such  as  instantaneously  switching,  zero  or 
infinite  resistance  transistor  switches  and  constant  air  gap 
flux.  Each  nodel  iaproved  a  particular  area  under  develop¬ 
ment  until  it  was  thought  to  be  a  suitable  representation  of 
an  actual  motor  or  power  conditioner  feature. 

The  first  model  incorporates  a  sguare  wave  of  back  emf 
(constant  flux)  and  endeavors  to  demonstrate  transistor 
commutation  effects.  The  transistor  switching  is  instanta¬ 
neous  and  is  accomplished  within  the  "commutation  and  Hall 
effect"  procedure  (see  Appendix  C)  . 

The  second  model  introduces  changes  in  the  transistor 
switching.  While  switching  is  still  instantaneous,  typical 
values  of  transistor  saturation  and  cutoff  resistance  are 
introduced,  and  commutation  is  controlled  by  the  commutation 
procedure  instead  of  being  accomplished  by  that  procedure. 
This  feature  leads  to  the  next  model  which  simulates  the 
actual  rise  time  and  saturation  delay  of  power  transistor 
switching.  This  model  is  the  basis  for  all  subsequent  simu¬ 
lations  involving  increased  model  complexities  such  as  the 
sinusoidal  flux  model. 

The  remaining  simulations  develop  the  voltage  and 
current  relationships  between  the  aotor  and  electronic 
commutator  final  drive  transistors  and  incorporate  sinu¬ 
soidal  aagnetic  flux  models  consisting  of  a  fundamental  and 
fifth  haracnic  as  described  in  Chapter  Four.  These  program 
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minus  IAB  equals  IB  and  ICA  ainus  IBC  equals  IC.  The  current 
flow  cannot  always  be  positive  as  depicted*  since  the  sub  of 
the  three  phase  currents  nust  be  identically  zero  by 
Kirchcff's  current  law.  Therefore*  at  least  one  of  the 
currents  IA*  IB  or  IC  aust  be  negative  at  any  given  tiae* 
representing  a  flew  of  current  returning  to  the  power 
supply.  From  the  current  wavefora  it  is  iaaediately 
apparent  which  phase  is  connected  to  the  negative  supply  and 
which  is  connected  to  the  positive  supply  siaply  by 
observing  the  polarity  of  the  current  (see  Figure  5.2).  The 
concept  of  the  positive  and  negative  supply  voltages  is 
totally  eguivalent  tc  a  single  positive  (or  negative)  supply 
as  used  in  [Bef.  1]  and  since  the  drive  currents  are  calcu¬ 
lated  froa  the  difference  of  two  fictitious  currents*  a 
single  supply  will  werk  in  the  siaulation*  but  the  split 
supply  enhances  the  syaaetry  of  the  notor  currents  and  volt¬ 
ages.  This  change  allows  for  iiproved  modeling  and  conceptu¬ 
alization  of  transient  studies  as  discussed  in  the  next 
secticn. 

C.  HCIEI  DEVELOP!! Ell 

The  basic  siaulation  of  the  brushless  DC  motor  by  Thomas 
used  a  single  positive  or  negative  power  supply  and  superim¬ 
posed  armature  leg  currents  to  produce  the  measured  motor 
performance.  Motor  drive  was  then  realized  by  multiplying 
averaged  armature  currents  by  a  torque  factor  for  the  model 
under  development.  In  the  this  model  the  supply  is  a  split 
supply  of  equal  voltages.  Armature  current  is  assigned  a 
positive  sign  if  it  flows  in  the  positive  direction  (i.e. 
into  the  motor)  as  shown  in  Figure  2.2  and  a  negative  sign 
if  it  flows  in  the  opposite  direction  (i.e.  out  of  the 
motor)  as  discussed  previously.  A  secondary  result  of  the 
model  of  commutating  transistors  and  split  supply  is  that 
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Bearranging  Eg  nation  2. 1  for  the  highest  derivative  on  the 
left  side  yields: 


dl/dt  *  (1/L)  Es  (t)  -  (B/L)I  (t)  -  (1/L)  KbWm  (t) 


(2.3) 


The  mechanical  torgue  balance  equation  is: 


Ts(t)  =  JdWa  (t)/dt  ♦  BWa(t)  ♦  Tl(t) 


(2.4) 


Eearranging  yields: 

dWa(t)/dt  *  (1/J)  Ta  (t)  -  (1/J)  Tl  (t)  -  (B/J)  Ha  (t)  (2.5) 

The  total  viscous  friction  (B)  and  total  load  inertia  (J)  is 
represented  as  follows: 


B  =  Ba  ♦  Bl 


(2.6) 


J  *  Ja  ♦  J1  (2.7) 

where  Bl  =  Blp/N  and  J1  *  Jlp/N.  Blp  is  the  viscous 

friction  associated  with  the  load*  Jlp  is  the  load  inertia 
and  N  is  the  gear  reduction  ratio.  The  link  between  the 
electrical  and  aechanical  eguations  is  the  torgue  constant* 
which  is  derived  froa  the  magnetic  flux: 


Tm  (t)  =  Ka$I  (t)  =  Ktl  (t)  . 


(2.8) 


The  basic  eguations  have  been  enhanced  in  this  xodel  to 
take  advantage  of  inherent  sign  conventions,  particularly  in 
direction  of  current  flow.  This  in  turn  Bakes  the  coaauta- 
tion  procedure  more  apparent  in  the  generated  graphics.  For 
example*  the  primary  aotor  currents  are  depicted  in  Figure 
2.2*  each  having  an  assigned  positive  direction  as 
pictured.  The  difference  between  IAB  and  ICA  yields  IA*  IBC 
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TABLE  II 

Sensor 

and 

Switching  Logic 

Counterclockwise 

Botation 

BCTOH 

BPS 

BPS 

BPS 

PHASE 

PHASE 

PHASE 

PCSITION 

i 

B 

C 

A 

B 

C 

0 

0 

1 

0 

open 

neg 

pcs 

30 

1 

1 

0 

pos 

neg 

open 

60 

1 

0 

0 

pos 

open 

neg 

90 

1 

0 

1 

open 

pos 

neg 

120 

0 

0 

1 

neg 

pos 

open 

150 

0 

1 

1 

neg 

open 

pos 

Clockwise  Notation 

BOTCB 

BPS 

BPS 

BPS 

PHASE 

PHASE 

PHASE 

POSITION 

A 

E 

c 

1 

B 

C 

0 

0 

1 

0 

open 

pos 

neg 

30 

1 

1 

0 

neg 

pos 

open 

60 

1 

0 

0 

neg 

open 

pos 

90 

1 

0 

1 

open 

neg 

pos 

120 

0 

0 

1 

pos 

neg 

open 

150 

0 

1 

1 

pos 

open 

neg 

C.  BC10B  SISTEH  EQUATIONS 

Derivation  of  tie  basic  notor  system  equations  (egua 
tions  2.1  through  2.8)  was  accomplished  by  Thomas  i 
[Bef.  1]  and  they  are  repeated  here  for  convenience.  Fro 
Kirchoff*s  voltage  lav,  the  stator  voltage  must  be: 

Es  (t)  =  Ldl/dt  ♦  BI(t)  ♦  Eb(t)  (2.1) 

where  the  back  emf  voltage  can  be  described  by: 


Eb(t)  =  (Km  $)Wm(t)  =  KbWm(t) 
or  Eb  (t)  =  Kbdfim  (t)  /dt 
since  d€m(t)/dt  =  Wm(t). 


(2.2) 
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Figure  2.3  (lotcr  Stator  and  Binding  Configuration 

rotation  is  tabulated  in  Table  II,  and  indicates  the  sensor 
logic  level  (RPS) ,  and  connection  of  the  phase  terninals  for 
both  clockwise  and  counterclockwise  direction  for  each 
thirty  degrees  of  rotor  position.  Note  that  regardless  of 
direction  of  rotation,  the  rotor  position  logic  remains 
constant  for  a  given  rotor  position,  and  that  rotor  posi¬ 
tions  between  180  degrees  and  360  degrees  are  indistingui¬ 
shable  from  rotor  positions  between  zero  and  180  degrees 
because  of  the  rotor  symmetry. 


Figure  2.2  Bctor  and  Connotating  Transistors 


Figure  4.4.  Since  tie  phase  windings  are  separated  by  sixty 
mechanical  degrees,  the  flux  in  each  phase  is  seen  to  have  a 
sixty  degree  phase  shift  with  respect  to  the  adjacent 
phases.  The  explanation  of  the  shape  of  the  model  flux  is 
discussed  in  Chapter  Four.  The  electronic  commutation  must 
he  performed  in  such  a  manner  as  to  complete  two  electrical 
cycles  fcr  each  mechanical  rotation  since  the  flux  distribu¬ 
tion  completes  two  cycles  in  one  rotor  rotation.  The  commu¬ 
tation  cf  the  armature  (stator)  current  is  accomplished 
every  thirty  degrees  cf  rotor  rotation  because  a  new  winding 
"enters"  the  magnetic  field  with  every  thirty  degrees  of 
rotation  simultaneously  with  another  winding  "leaving"  the 
magnetic  field.  The  sensor  and  switching  logic  for  motor 


TABLE  I 

Typical  Botor  Paraaeters 


Parameter 

Stater  resistance,  Ra 

Stator  Inductance,  la 

Torque  coefficient,  Kt 

Reverse  voltage 

coefficient,  Kb 

Ector  inertia,  Jm 

Viscous  fpiction 

coefficient,  Em 


1A ias 

Onits 

1.37 

Ohms 

0.0016 

Henrys 

15.9 

Oz-In/amp 

0.112 

Volt-Sec/Rad 

0.001 

Oz-In/Sec 

0.022 

Oz-In-Sec/Rad 

B.  HCTOE  COHFIGOH ATICB 

The  motor  simulated  in  this  thesis  is  a  three  phase, 
four  pole  device  typical  of  modern  brushless  DC  motors  as 
described  in  [Bef.  1],  Typical  values  for  the  motor  parame¬ 
ters  are  contained  in  Table  I.  The  listed  values  of  resis¬ 
tance  and  inductance  are  as  measured  across  two  windings  as 
the  motor  is  only  a  three  wire  device  as  pictured  in  Figure 
2.2.  It  should  be  noted  that  the  motor  inertia  is  very  low; 
a  typical  rotor  diameter  is  of  the  order  of  one  inch  and 
results  in  a  very  small  mechanical  time  constant.  A  tran¬ 
sverse  cross  sectional  view  of  rotor  and  stator  for  a  the 
model  motor  is  shown  in  Figure  2.3,  with  the  rotor  in  the 
zero  (or  180)  degree  position.  The  magnetic  rotor  can  be 
seen  to  have  north  and  south  poles  separated  by  ninety 
degrees.  This  generates  a  flux  distribution  which  alternates 
polarity  every  ninety  mechanical  degrees  as  depicted  in 


saturation  and  finally  drive  it  off  quickly  using  reverse 
drive.  This  can  he  accomplished  neatly  and  siaply  by 
bypassing  the  base  resistor  with  a  suitably  sized  capacitor 
as  discussed  in  [Bef.  6  :p.  784]. 


C.  BCD El  OF  SNITCHING  TRANSISTOR  DYNAMICS 

The  preceding  discussion  leads  into  the  modeling  of  the 
transistor  switches  as  dynamic  elements.  The  most  desirable 
transistor  model  would  be  of  low  order  and  not  add  signifi¬ 
cantly  tc  the  program  complexity  since  the  motor  dynamics 
are  the  prime  concern. 

The  initial  models  were  developed  with  the  transistor 
switches  acting  instantaneously  between  cutoff  resistance 
and  saturation  resistance.  This  was  satisfactory  as  long  as 
the  transients  generated  remained  small.  Osing  the  variable 
step  size  integration  routines  available  with  CSMP,  it  was 
found  that  the  rapid  switching  drove  the  integration  step 
size  minutely  small  ard  computer  time  for  the  simulation  to 
reach  steady  state  became  excessively  long.  Based  upcn  the 
discussion  given  in  the  preceeding  section  on  transistor 
switching,  an  exponential  rise  and  fall  was  developed,  and 
both  overdrive  and  reverse  drive  simulated  by  clipping  the 
exponential  rise  and  decay  within  limits  of  the  saturation 
resistance  and  cutoff  resistance.  This  not  only  provided  a 
good  simulation  of  the  transistor  switching  with  user  vari¬ 
able  parameters,  but  also  eased  the  computer  time  problem 
somewhat.  The  resulting  switching  characteristics  are  shown 
in  Figure  3.3,  (with  exaggerated  tine  constants  for  clarity) 
and  ccmpare  favorably  with  those  shown  in  Figure  3.2  (a)  and 
(b)  .  Figure  2.2  suggests  that  the  transistors  are  bipolar 
NPN  devices.  This  was  done  for  pictorial  reference  only. 


TABIE  III 

Son*  Power  Transistor  Specifications 


Device 

Ic  (flax) 
&aB§ 

B  (sat) 
Ohms 

ts 

s 

Diss 

Watts 

Vbr  (sus) 
Volts 

M1M45N15 

(HCS) 

45* 

.061 

250 

150 

MTM40N18 

(MOS) 

40* 

.08% 

250 

180 

MJ16016 

(Bipolar) 

20 

2.2 

250 

450 

M  J 1501 1 
(Bipolar) 

10 

200 

250 

2H3430 

(Bipolar) 

5 

.2 

4.0 

150 

100 

2N2777 

(Eipclar) 

30 

.06 

14** 

200 

150 

2N6250 

(Eipclar) 

15 

3.5 

175 

275 

2B6339 

(Eipclar) 

25 

1.0 

200 

120 

2S1820 

(Eipclar) 

15 

.1 

25** 

250 

250 

S1C2231 

(Eipclar) 

30 

.08 

6* 

200 

200 

♦  ID 

SB  (on) 

**ts+tf 

#td+tr 

These  switching  transistors  may  he  a  pair  of  complementary 
symmetry  Class  B  amplifier  output  devices  to  take  advantage 
of  emitter-follower  characteristics  or  might  be  MOS  devices 
to  realize  high  transistor  power  efficiency.  The  model  as 
presented  is  accurate  for  bipolar  transistor  dynamics  hut 
because  switching  time  constant  and  on  and  off  resistance 
are  the  cnly  parameters  for  the  model,  it  should  be  exten¬ 
dable  to  include  MOS  devices  as  well.  The  parameters  used 
for  the  simulation  program  were  a  switching  time  constant  of 


1.5  Microseconds  and  a  saturation  resistance  of  .05  ohms. 
These  values  were. chosen  as  more  or  less  ideal  values,  the 
combination  of  which  would  probably  not  be  realized  in  any 
one  device.  In  Table  III  some  transistor  specifications  are 
listed  for  devices  which  might  be  desirable  to  use  for  this 
motor  controller  application.  The  fast  switching  times  for 
the  bipolar  transistors  are  their  salient  feature,  while  the 
low  "on”  resistance  of  the  HOS  devices  makes  them  desirable. 


TRANSISTOR  SATURATION  MODEL 


ROTATION  ANGLE,  DEGREES 

Figure  3.3  Transistor  Model  Saturation  Characteristics 

The  model  developed  to  this  stage  has  ignored  the  effect 
of  peak  reverse  voltage  across  the  transistors  that  occurs 
at  the  moment  of  commutation  due  to  the  opening  of  the 
inductive  circuit.  To  provide  the  required  protection. 


damper  diodes  are  used  to  clasp  the  spake  voltage  tc  the 
positive  or  negative  supply.  The  diodes  not  only  provide 
overvoltage  protection  for  the  transistors,  but  oust  handle 
nearly  the  full  value  of  steady  state  current  for  a  brief 
period  of  tine  as  well.  The  effect  of  the  diode  clanping 
can  be  seen  innediately  by  coaparing  the  figures  depicting 
notor  current  with  and  without  the  diodes  in  the  simulation. 
(See  Figures  5.1  and  5. 2)  Note  that  the  current  dip  is  mark¬ 
edly  reduced  when  the  diode  is  in  the  circuit.  The  diodes 
provide  a  commutation  of  the  transistor  current  and  maintain 
current  flow  in  the  notor  during  the  switching  of  the  drive 
current  between  two  phases.  Additional  information 
concerning  diode  comcutation  was  found  in  £Be£.  9  :p.  371] 
concerning  the  so-called  "Jones"  commutation  circuit  for  a 
EC  meter  controller. 

The  diode  model  is  very  simple  and  features  instanta¬ 
neous  switching  between  saturation  and  cutoff  resistance 
identical  in  value  to  that  of  the  transistors.  This 
approach  did  not  add  any  appreciable  complexity  to  the  inte¬ 
gration  problem  because  the  diodes  pick  up  the  current  as 
the  transistor  turns  off,  smoothing  out  the  transient 
current  instead  of  generating  another  current  transient. 
The  switch  of  the  current  to  the  diode  is  triggered  by  the 
voltage  spike  caused  by  the  transistor  turning  off,  similar 
to  what  would  be  expected  in  the  actual  circuit.  The  lack  of 
a  diode  voltage  drop  before  the  knee  of  the  diode  conduction 
curve  is  the  only  addition  that  might  be  made  to  the  model, 
but  was  deemed  unnecessary  as  the  voltage  drop  across  the 
diode  at  full  meter  current  would  still  be  less  than  one 
volt  and  of  the  same  erder  of  magnitude. 


IV.  MAGNETIC  FIOI  MODEL 


A.  SIHMETBIC  FIFTH  BABHONIC  FIOX 

Thomas  [Bef.  1]  derived  a  combination  fundamental  and 
fifth  harmonic  reverse  voltage  waveform  from  empirical  meas¬ 
urements.  This  model  has  been  used  in  this  thesis  with 
modifications.  The  phase  relationship  of  the  fifth  harmonic 
to  the  fundamental  has  been  altered  to  provide  the  same 
shape  voltage  waveform  across  two  windings  instead  of  across 
only  one.  This  agrees  with  the  empirical  data  and  also  gives 
rise  to  a  new  magnetic  flux  model,  the  shape  of  which  agrees 
in  form  with  traditicral  DC  machines  [Ref.  8].  The  general 
shape  of  the  flux-density  distribution  of  a  DC  motor  is 
shown  in  Figure  4.1. 


co^A  I 


Figure  4.1  Flux  Density  Distribution  due  to  Botor  Field 


The  dip  in  the  center  of  the  flux-density  distribution  can 
test  fce  explained  by  considering  Figures  4.2  and  4.3  which 
are  flux  maps  of  the  magnetic  field  due  to  the  rotcr  magnets 
and  stater  current  respectively  and  confirm  the  distribution 
shown  in  Figure  4.1  for  either  stator  or  rotor  considered. 
In  Figure  4.2  and  Figure  4. 3,  separation  of  flux  lines  indi¬ 
cates  degree  of  field  strength,  with  greater  separation 
indicating  less  field  strength. 
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Figure  4.2  Flux  Bap  with  only  Rotor  Field  Considered 
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Figure  4.3  Flux  Bap  with  only  Stator  Field  Considered 
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The  model  flax-density  distribution  is  depicted  in 
Figure  4.4  showing  the  phase  relationship  between  the  three 
windings.  In  Figure  4.5  the  phase  'A*  flux  is  shown  along 
with  the  resulting  generated  voltage  between  terminals  *A' 
and  *Bf.  The  dotted  lines  in  each  figure  represent  the 
commutation  periods  during  which  each  winding  is  connected 
to  the  positive  supply,  the  negative  supply,  or  effectively 
disconnected  from  supply  by  the  switching  power  transistors. 
The  ccmmutation  effects  shown  by  the  dotted  lines  in  these 
figures  represent  motor  currents  for  counterclockwise  (posi¬ 
tive)  rotation  as  the  product  of  the  flux  and  the  motor 
current  is  seen  to  always  be  positive.  Since  the  motor 
torque  is  derived  from  the  product  of  the  magnetic  flux 
density  distributions  and  the  current  flowing  in  each  phase 
winding,  this  results  in  a  positive  motor  torque  and  rota¬ 
tion  in  the  positive  direction.  The  composite  flux  density 
distribution  for  the  rotating  motor  is  displayed  in  Figure 
4.6,  and  is  used  to  derive  the  normalizing  constant  for  the 
flux  density  by  maintaining  the  average  value  for  the  torque 
coefficient  at  the  constant  value  for  the  simple  models. 
This  composite  flux  density  distribution  is  seen  to  have  an 
average  value  equal  to  the  torgue  constant  of  the  simple 
models,  when  the  conversion  from  webers  to  inch-ounces  per 
amp  is  applied: 

(in-cz/amp)  x  (.0C7Q7)  =  webers  (4.1) 

This  figure  also  provides  a  quick  check  on  the  back  emf 
constant  since  the  numerical  value  of  both  constants  must  be 
identical  in  WKS  units. 

Tc  demonstrate  clockwise  rotation,  the  flux  model 
remains  fixed  with  respect  to  the  shaft  position,  but  the 
sign  cf  the  currents  must  be  reversed.  This  amounts  to 
reversing  the  polarity  of  the  commutation  scheme  for  a 
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Figure  *1.4  Hodel  Flu* — Phase  Belationship 
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Figure  4.5  Bodel  Flu*  and  Back  EHF  Voltage 
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MAGNETIC  FLUX  MODEL 

COMPOSITE  FLUX  VARIATION 


Figure  4.6  Model  Composite  Flux 

counterclockwise  rotating  motor.  As  can  be  seen  from  Table 
II  this  is  the  same  sequence  of  commutation  as  for  counter¬ 
clockwise  rotation  advanced  by  ninety  degrees,  which  should 
be  expected  from  the  rotor  geometry  since  the  south  pole  of 
the  rotor  is  displaced  ninety  degrees  from  the  north  pcle  of 
the  rotor. 


E.  FIFTH  HARMONIC  FIUI  HIT H  DRIVE  CURRENT  INTERACTION 

The  interaction  of  the  two  magnetic  fields,  one  due  to 
the  rotor  magnets  and  the  other  due  to  the  stator  current 
causes  distortion  of  the  air  gap  flux  and  hence  the 
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generated  voltage  waveform.  Eeference  eight  indicates  that 
such  an  interaction  is  expected,  but  as  to  what  degree  the 
flux  density  distribution  and  the  reverse  voltage  waveform 
would  be  distorted  is  unknown.  The  expected  effect  of  the 
two  magnetic  fields  interacting  is  to  skew  the  symmetrical 
distribution  as  is  shewn  in  Figure  4.7  It  is  realixed  that 
the  profile  of  the  flux  model  of  Figure  4.5  and  the  total 
flux  density  distribution  as  shown  in  Figure  4.7  do  not 
agree,  however  the  flux  model  shows  a  high  degree  of  corre¬ 
lation  to  empirically  measured  experimental  test  data  for 
this  type  of  motor  £  Bef .  11]. 


Figure  4.7  Flux  Density  Distribution 
Demonstrating  Interaction. 

The  skewed  shape  cf  the  resultant  flux-density  distribu¬ 
tion  due  to  stator  current  interaction  can  be  simulated  by 
slipping  the  phase  relationship  of  the  fifth  harmonic  with 
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respect  to  the  f undamental.  If  the  phase  shift  is  singly  a 
constant,  then  two  windings  will  produce  essentially  the 
sane  shape  of  back  IB?  voltage  as  the  previous  model,  but 
combination  with  the  other  winding  will  not  be  similar  or 
symmetric.  Since  current,  as  modeled,  essentially  flows  in 
two  windings  at  a  tine,  the  distortion  must  be  a  function  of 
the  current  and  the  third  winding  flux  must  remain  undis¬ 
torted  as  long  as  that  winding  has  no  current  flowing. 
This  thinking  led  to  the  model  for  interactive  flux  depicted 
in  Figure  4.8.  The  phase  relationship  between  fundamental 
and  harmonic  is  identical  to  the  non-interactive  model  if 
the  drive  current  is  zero.  In  the  presence  of  drive 
current,  the  fifth  harmonic  is  shifted  smoothly  and  continu¬ 
ously  until  at  maximum  current,  the  phase  shift  is  thirty 
electrical  degrees.  Even  with  this  visible  warping  cf  the 
flux,  the  back  emf  voltage  waveform  is  recognizeatle  as  a 
fundamental  and  a  fifth  harmonic,  and  this  effect  might  be 
smaller  than  depicted,  and  of  little  consequence  to  the 
motor  designer.  Additionally,  no  model  of  magnetic  satura¬ 
tion  has  been  implemented,  which  would  have  the  effect  of 
flattening  the  flux  peaks  and  giving  a  more  constant  maximum 
flux  and  hence  torque,  resulting  in  smoother  motoring.  Such 
effects  are  expected  in  an  actual  motor,  but  at  what  value 
of  current  the  saturation  should  begin  will  have  to  be 
determined  by  measurement  of  a  particular  motor. 
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Figure  4.8  Hodel  Flux  and  Back  EHF  Vgltage 
fieaonstrating  Drive  Current  Interaction 
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With  minor  changes  to  the  program,  different  magnetic 
flux  shapes  can  be  simulated  and  motor  performance  can  be 
compared.  For  example,  if  the  stator  was  designed  with 
interpoles,  the  flux  would  take  on  a  flatter  shape  across 
the  tcp  of  the  flux  curve,  which  would  in  turn  give 
different  characteristics  to  the  back  emf  voltage,  which 
could  be  empirically  measured.  If  the  torque  coefficient  of 
the  motor  was  known,  (it  can  be  measured)  the  flux  distribu¬ 
tion  could  be  deduced  as  it  has  been  for  these  simulations 
and  the  complete  motor  performance  curves  predicted.  The 
simulations  give  a  good  representation  of  transient  and 
steady  state  motor  and  transistor  currents  [Ref.  10],  from 
which  voltage  and  current  stresses  applied  to  the  transistor 
junctions  can  be  determined.  The  results  of  a  simulation 
for  a  desired  steady  state  load  may  be  analyzed  to  size  the 
transistors  and  diodes  required  for  a  specific  application. 
Given  the  heat  transfer  characteristics  of  a  motor  and/or 
heat  sink  assembly,  the  thermal  rise  of  the  power  condi¬ 
tioner  and  the  motor  stator  may  be  predicted  by  implementing 
a  simple  heat  transfer  model. 

From  the  literature  [Hef.  3]  it  was  noted  that  when  some 
thirty  degrees  of  commutation  advance  was  used  to  obtain 
maximum  performance  from  the  motors  being  tested,  an 
increase  of  2307.  over  rated  horsepower  was  realized  for 
short  periods.  This  simulation  program  has  the  capability 
of  shifting  the  commutation  angles  to  realize  extended 
performance.  This  is  left  for  future  development  of  motor 
performance  studies. 
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CENTER  POINT  TRAJECTORY 


o 


ROTATION  ANGLE,  DEGREES 


Figure  5.8  Trajectory  of  Center  Connection 

the  variation  in  speed  would  be  reduce!  to  the  order  of  five 
percent  due  to  tfce  increased  tine  constant,  (J/E) . 
Additionally,  the  speed  ripple  decreases  significantly  if 
the  applied  voltage  is  increased  with  no  increase  in  load. 
Uith  a  load  toriue  of  64  inch  ounces  and  a  power  supply 
voltage  of  sixty  volts,  speed  ripple  is  hardly  noticeable  at 
the  steady  state  speed  of  approximately  4800  rpm  due  tc  the 
increased  angular  momentum. 
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resulting  in  smoother  rotation.  The  design  tradeoff  to  this 
would  he  that  only  half  of  the  supply  voltage  would  be 
available  across  any  one  winding  with  the  center  grounded. 
If  the  center  terminal  is  left  to  float,  the  voltage  across 
two  windings  is  the  full  supply  voltage,  but  the  voltage 
applied  across  a  single  winding  and  the  current  through  that 
winding  is  a  functicn  of  the  back  eaf  voltage  generated  in 
both  windings.  Since  the  back  emf  voltage  is  varying  sinu¬ 
soidally  in  each  winding  and  one  winding  is  120  electrical 
degrees  cut  of  phase  with  respect  to  each  of  the  other  two 
windings,  a  larger  proportion  of  the  supply  voltage  may  be 
applied  to  one  of  the  conducting  coils  than  the  other.  The 
current  through  the  cne  coil,  being  equal  to  the  sum  of  the 
currents  in  the  other  two  coils,  is  at  a  higher  value  when 
the  air  gap  flux  is  greatest  than  if  the  center  connection 
were  grounded,  and  will  provide  more  motor  torque.  The 
trajectory  of  the  center  connection  is  shown  in  Figure  5.8 
and  can  be  compared  tc  the  single  phase  flux  characteristic 
in  Figure  4.4  and  the  composite  flux  variation  in  Figure  4.6 
to  see  that  the  voltage  aids  the  current  flow  in  each 
winding  as  it  cones  under  the  peak  of  the  magnetic  flux,  as 
the  voltage  at  the  center  connection  is  negative  at  during 
the  positive  flux  peaks,  and  positive  during  the  negative 
flux  peaks,  aiding  current  flow  when  compared  to  a  grounded 
center  terminal.  The  relative  scale  on  the  ordinate  of  this 
figure  indicates  that  the  effect  is  proportional  tc  motor 
speed. 

The  fluctuation  of  motor  speed  demonstrated  ty  this 
model  has  been  of  the  order  of  ten  percent  of  peak  motor 
speed,  somewhat  more  than  in  Thomas*  simulations  [Bef.  1]. 
This  can  be  attributed  to  the  combination  of  the  transistor 
switching  time  constants,  the  current  rise  time  constant 
(1/R)  and  extremely  small  rotor  inertia.  If  the  load 
inertia  was  of  the  order  of  magnitude  of  the  rotor  inertia. 
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10  these  figures,  a  demand  for  reversal  is  initiated  at 
.02  seconds.  Since  tie  rotational  inertia  of  the  rotor  keeps 
it  rotating  in  the  positive  direction  for  a  short  time  after 
reverse  current  begins  to  flow,  a  commutation  point  is 
reached  and  the  reversal  current  is  switched  between  phases 
A  and  C.  As  the  rotor  stops  and  the  motor  reverses  under  the 
influence  of  negative  torque,  the  commutation  point  is  again 
reached  tut  from  the  opposite  direction  and  the  current  is 
switched  from  phase  C  to  A.  There  is  a  large  current  tran¬ 
sient  in  the  figure  at  this  point  which  is  also  seen  in 
figure  5.7  as  a  transient  positive  torque,  probably  due  to  a 
combination  of  the  instantaneous  nature  of  the  diode  model 
and  the  transistor  time  constants.  The  reality  of  such  a 
current  and  torque  spike  is  questionable,  but  if  the  dynamic 
nature  of  the  transistors  and  diodes  used  in  the  power 
conditioner  are  not  similar,  such  undesirable  transients  may 
be  generated,  even  to  the  point  of  punch  through  of  a  junc¬ 
tion.  The  effect  of  the  torque  spike  is  seen  in  Figure  5.6 
as  a  slight  interruption  of  the  otherwise  smooth  transition 
between  forward  rotation  and  reverse  rotation. 

The  simulation  can  be  a  powerful  analytic  and  design 
tool.  For  example,  the  addition  of  damper  diodes  was  seen 
to  maintain  current  flow  during  transistor  switching  (demon¬ 
strated  in  Figures  5.1  and  5.2)  which  has  the  effect  of 
maintaining  a  higher  average  torque  with  diodes  than  without 
as  has  been  previously  discussed.  It  follows  that  the 
faster  the  diodes  operate,  the  better  the  motor  performance 
which  can  be  realized  in  this  context. 

The  availability  of  the  center  of  the  three  windings  for 
use  as  a  ground  terminal  may  prove  to  be  advantageous  in 
some  applications.  The  current  return  path  to  the  power 
supply  provided  by  the  center  connection  would  not  be  opened 
by  the  switching  transistors,  and  the  current  flow  would  not 
be  interrupted  in  the  unswitched  winding  during  commutation. 
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Figure  5.4  Phase  B  Current  shoving  Reversing  Transient 


Motor  reversal  at  speed  can  be  simulated  under  varying 
load  conditions.  In  the  simulation  the  load  torque  opposes 
motor  rotation,  and  so  helps  the  motor  to  slow  to  a  stop, 
then  opposes  the  reversal.  If  a  viscous  load  was  simulated 
(as  in  driving  a  dynamometer),  slowing  to  a  stop  would  take 
longer  because  the  lead  would  decrease  with  speed,  and  the 
rise  to  steady  state  in  the  reverse  direction  would  be 
quicker  due  to  the  reduced  load  at  low  speed.  In  general, 
this  type  of  reversal  is  characterized  by  large  transient 
motor  currents  since  the  back  emf  voltage  is  initially  of 
polarity  that  aids  current  flow  instead  of  inhibiting  it. 
Ihe  reverse  torque  generated  by  this  current  will  quickly 
slow  the  motor  to  a  stop  where  there  is  no  back  emf  gener¬ 
ated,  and  as  the  rotation  reverses  the  generated  voltage 
changes  polarity  and  once  again  opposes  current  flew.  The 
reversal  effects  can  be  seen  in  the  current  waveforms  of 
Figures  5.3,  5.4,  5.5  and  5.6  which  depict  motor  reversal 
(from  counterclockwise  rotation  to  clockwise  rotation)  at 
time  -  . C2  seconds. 


Figure  5.2  Hotor  Current  with  Diode  Confutation 
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Pigure  5.1  ilotor  Current  without  Diode  Couuutation 


The  final  versions,  in  which  the  transistor  and  diode 
models  and  sinusoidal  back  emf  voltage  are  incorporated  are 
used  to  characterize  an  actual  aotor  and  power  conditioner. 
The  performance  curves  provided  by  Thomas  [Bef.  1]  were  used 
for  comparison.  As  a  result,  all  simulations  presented  here 
used  a  net  supply  voltage  of  thirty  volts,  or  plus  and  minus 
fifteen  volts  applied  to  the  transistors.  As  stated  earlier 
the  object  of  this  thesis  was  to  provide  a  tool  to  evaluate 
motor  and  switching  transistor  performance  for  various 
design  parameter  values  and  operating  conditions.  Seme 
deviation  from  these  performance  curves  was  anticipated,  hut 
overall  characteristics  of  the  torque,  current,  and  power  vs 
speed  curves  remain  of  the  same  form. 

At  the  time  of  this  writing,  provisions  are  being  made 
to  gather  mere  comprehensive  test  data  from  a  typical  metor. 
It  is  well  understoed  that  only  a  limited  number  of  the 
model  waveforms  may  be  confirmed  by  measurement.  The  brush¬ 
less  DC  motor  is  a  three  wire  device,  and  the  voltage  across 
any  twe  windings  may  be  conveniently  measured.  The  simula¬ 
tion,  however,  calculates  and  graphs  not  only  the  voltage 
across  each  pair  of  windings  but  in  addition,  the  voltage 
across  each  single  winding.  The  individual  currents  IA,  IB, 
and  IC  nay  also  be  measured  with,  for  example,  an  inductive 
coupler  or  by  measuring  the  voltage  drop  across  a  resistor. 
These  measured  curves  can  be  readily  compared  with  model 
results  and  model  parameters  adjusted  to  bring  about  gcod 
agreement. 


f.  HOPE!  PEBPORflANCE 


The  performance  cf  the  first  three  models  is  net  inter¬ 
esting  to  the  motor  designer  as  the  simplifications  involved 
in  these  simulations  reduce  the  resulting  behavior  to  text¬ 
book  PC  motor  dynamics.  As  a  development  tool  for  other 
simulations,  the  structure  may  be  of  some  interest.  These 
simulations  are  included  in  Appendix  C  for  this  reason.  The 
transistor  switching  of  the  fourth  and  all  subsequent  models 
is  of  interest  in  that  the  transistor  impedance  affects  the 
(L/R)  time  constants  of  the  motor  current  and  is  an  impor¬ 
tant  feature  in  understanding  the  behavior  of  the  simula¬ 
tions.  Because  of  the  varying  time  constants,  decay  of 
current  in  one  phase  and  buildup  of  current  in  another  phase 
at  commutation  is  net  compensatory.  The  current  decays 
faster  in  the  turn-off  winding  than  the  corresponding 
buildup  in  the  turn-cn  winding,  while  the  current  in  the 
third  winding  is  ostensibly  constant.  The  third  winding 
current  cannot  be  constant  because  all  three  currents  must 
sum  to  identically  zero.  The  result  is  the  sharp  dip  in  the 
center  of  the  current  waveform  of  the  unswitched  winding. 
Excellent  correlation  of  this  effect  is  seen  in  the  experi¬ 
mental  data  as  shown  in  Figure  (12)  of  [Ref.  10]  and 
[Ref.  11].  In  simulations  which  include  a  diode  model,  the 
dynamic  equations  that  determine  the  (L/E)  time  constant  are 
expanded  to  include  the  diode  model.  The  inclusion  cf  the 
diode  impedance  in  the  dynamic  time  constant  results  in  a 
"bracing  up"  of  the  current  waveform  of  the  unswitebed 
winding  yeilding  a  higher  average  motor  current.  This  can  be 
observed  by  comparing  Figures  5.1  and  5.2.  Models  with 
simulated  diodes  demonstrate  somewhat  more  torque  and  speed 
with  slightly  less  speed  ripple  than  those  without  diodes. 
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VI.  C0HCL05I0B  AMD  BECOMBBHDATIOBS 


A.  S09BAEI  AHD  COHCLOSIOH 

The  development  cf  a  detailed  transient  model  for  simu¬ 
lation  cf  the  instantaneous  performance  of  a  small  four  pole 
brushless  DC  motor  supplied  by  a  transistorized  power  condi¬ 
tioner  has  been  presented.  This  model  was  part  cf  a  compre¬ 
hensive  development  cf  a  fin  control  actuating  system  model 
for  the  simulation  cf  the  dynamic  performance  cf  a  cruise 
missile  fin  controller. 

The  results  of  the  numerical  simulation  of  the  instanta¬ 
neous  meter  and  power  conditioner  voltages  and  currents 
demonstrated  a  great  deal  of  correlation  to  reported  experi¬ 
mental  test  data  during  steady  state  operation.  The  dynamic 
reversal  of  the  model  was  performed.  The  corresponding 
dynamic  currents  predicted  by  the  model  were  analyzed  and 
large  transients  were  demonstrated  as  would  be  expected  in 
an  actual  motor.  The  conclusion  that  the  model  is  a  good 
representation  of  an  actual  motor  and  power  conditioner 
drive  transistors  is  supported  by  all  currently  available 
motor  data  and  analytical  extensions  thereof. 

B.  BECCHHEHDATIOHS 

This  model  will  he  further  improved  by  the  addition  of 
delay  between  the  Hall  sensors  and  the  commutation  algorithm 
to  more  realistically  simulate  the  electronic  logic  cf  the 
power  conditioner.  A  routine  for  the  calculation  cf  tran¬ 
sistor  and  diode  power  dissipation  and  possibly  motor  stator 
heating  will  be  included.  This  will  allow  maximum  perform¬ 
ance  to  be  extracted  from  the  motor  and  power  conditioner 
for  short  periods  of  time  without  excessive  heating. 


Additional  work  cn  the  dynaaic  aagnetic  flux  profile, 
particularly  in  the  area  of  aagnetic  saturation  of  the 
stator  iron  needs  tc  be  investigated  and  iapleaented  to 
accurately  siaulate  high  power  operation. 


SPACE  DEPENDENT  FLOX  MODEL 


Because  of  the  length  and  complexity  of  the  CSMP  model, 
a  simpler  model  was  designed  using  Digital  Simulation 
language  (DSL) .  Tie  advantages  to  DSL  include  complete 
compatibility  with  CSMP  and  a  more  or  less  interactive 
running  capability  which  the  NPS  installation  of  CSMP  lacks. 
If  the  complete  model  had  been  written  in  DSL,  the  advantage 
of  guick  turnaround  of  the  simulation  would  have  been  lest 
due  to  the  running  tine  of  the  program.  Ihe  NPS  version  of 
DSL  dees  not  have  tie  STIFF  integration  routine  used  with 
CSMP  and  the  Bunge-Kutta  integration  routine  in  both 
languages  had  some  computation  problems  with  the  fast  time 
constants  of  the  power  transistors.  The  routines  contained 
in  this  appendix  were  developed  using  [Ref.  14  and  15]. 

TITLE  THESIS  PROJECT 

TITLE  BBCSH1ESS  DC  HCTCR  POWER  CONDITIONER 

STORAG  Z 1  (500)  ,  Z2  (500)  ,Z3  (500)  ,Z4  (500)  ,  Z5  (500)  ,  Z6  (500) 

CONSI  PI  =  3.14159,  VIF=15.t  VIB  =  -15. 

******************** ******************* 

DYNAMIC 

*************************************** 

THETA  =  RAMP  (0.01 

THITAD  *  THETA ♦180. /PI 

EEHFA  =  (3.  *S  IN  (2.*TH2TA*  (  1 1*PI/6)  ) 

1  +  .  59*SIN(1 0. *THETA+ (1  *PI/6)  ) ) 

EEMFE  =  (3.  *SIN  (2.  *THET  A+  (7*PI/6)) 

1+.59*SIN  (10.*TEETA+  (5*PI/6)  )) 

BEEFC  “  (3.  *SIN  (2. *THETA+  (3*PI/6)  ) 

1  +  .59*SIN ( 1 0. *THETA+ (9*PI/6)  )) 

AEXIFF=  BEMFA  -  EEMFB 

BCDIFF=  B2MFB  -  EEMFC 

CADIFF=  BEMFC  -  BEMFA 

IF  (THETA D  .GE.  180.)  GO  TC  45 
THCON  =  THETAD 
GO  TO  46 

45  THCCN  =  THETAD-  180. 

46  CONTINUE 

IF  (THCON  .GE.  0..AND.  THCON  .LT.  30.)  GO  TO  50 

IF  THCON  .GE.  30. .AND.  THCON  .LT.  60.'  GO  TO  51 

IF  THCON  .GE.  60. .AND.  THCON  . LT.  90.  GO  TO  52 

IF  (THCON  .GE.  90. .AND.  THCON  .LT.120.  GO  TC  53 

IF  THCON  . GE. 1 20. .AND.  THCON  .IT. 150.  GO  TO  54 


/6>) 


1  +  .59*SIN  (10.*TEETA+ (5*PI/6)  )) 
BEEFC  =(3.  *SIN  (2.  +  THETA+ (3*PI/6) 
1  +  .59*SIN  (1 0.  *THETAM9*PI/6)  )) 


THCON  .LT.  30. 
THCON  .LT.  60. 
THCON  .LT.  90. 
THCON  .LT.120. 
THCON  . LT. 1 50. 


IF  (THCON 

srf  =  o. 

SW2  =  0. 
SW3  =  0. 


GE. 150. .AND.  THCON  .LT.180.’ 


GO  TO  50 
GO  TO  51 
GO  TO  52 
GO  TO  53 
GO  TO  54 
GO  TO  55 


SH€  =  0. 

BEBFT  =  BEMFC  -  EEMFB 
VN 1  =  VI F  -  VEHFC 
VN2  =  VIB  ♦  VEBFB 
GO  TO  60 
SHI  =  1. 

SH2  =  0. 

SH3  =  0. 

SHU  =  1. 

SH5  =  0. 

SH6  =  0. 

BEBFT  =  BEHF A  -  EEMFB 
VN  1  =  VI F  -  VEMFA 
VN2  =  VIB  +  VEBEB 
GO  TO  60 
SHI  =  1. 

SH2  =  0. 

SH3  =  0. 

SHU  =  0. 

SH5  =  0. 

SH€  =  1. 

BEBFT  =  BEMFA  -  BEMFC 
VN  1  =  VI F  -  VEMFA 
VN2  =  VIB  ♦  VEMIC 
GO  TO  60 
SHI  =  0. 

SH2  =  0. 

SH3  =  1. 

SF4  =  0. 

SH5  =  0. 

S¥  6  =  1. 

BEBFT  =  BEMFB  -  EEMFC 
VN  1  =  VI F  -  VEMIB 
VN2  =  VIB  ♦  VEBFC 
GO  TO  60 
SHI  =  0. 

SH2  =  1. 

S13  =  1. 

SH4  =  0. 

SH5  =  0. 

SH  €  =  0. 

BEBFT  =  BEMFB  -  BEMFA 
VN  1  =  VI F  -  VEMFB 
VN2  =  VIB  ♦  VEBFA 


VN2  =  VIB  ♦  VEBFA 
GO  TO  60 
55  SHI  =  0. 

SH2  =  1. 

SH3  =  0. 

SW4  =  0. 

SH5  =  1. 

SH€  =  0. 

BEBFT  =  BEMFC  -  BEMFA 
VN  1  =  VI F  -  VEMFC 
VN2  =  VIB  ♦  VEBFA 
60  CONTINUE 

IA  =  4.*  f SH1-SH2) 

IB  =  4.*  (SH3-SH4) 

IC  =  4.  *  (SH5-SH6) 

**************  ************************* 
SAMPLE 

*************************************** 

*  GET  TEE  NUMBEBS 


SW1-SH2) 

SH3-SH4) 

SH5-SH6) 

************************ 


THETAD 
BEMFA 
BEMFB 
BEMFC 
IA 


Z6  (I)  =  IB 
Z7jl)=  IC 
Z8jlj  =  BEMFT 
WHITE  (30,20 
WRITE  31,20 
WHITE  32,20 
WHITE  33,20 
WHITE  34,20 


WHITE  32,20' 
WHITE  33,20; 
WHITE  34,20 
WHITE  35,20 
WHITE  36,20 
WRITE  37-20, 
rOBHAT  F20.6 


ll\l\ 

z3  ;i 
Z4  I 
Z5  I 
Z6  I 
Z7  I 
Z8  I 


1=1*1 

COHTB1  FIHTIM=6.28,  EiIT=.01,  DELS=.01 

BHD 

STOP 
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Beam  DISSPLA  prograa: 

y  REAL  Z1  (629r,Z2  (629)  ,Z3  (629)  ,Z4  (629)  ,Z5  f629)xZ6  (629) 
BEAI  z7  (629)  ,Z£  (629)  ,  Z9  (629)  , Z  10  (629/ ,  Z 11  (629) 


1=0  . 

1  1=1+1 

BEAT  (30,5, END  =  3)  Z1  (I) 

BEAD  3lj5;END  =  3  Z2  I 
BEAD  32,5, END=3  Z3  I 
BEAD  33,5,END  =  3;  Z4  (I 
READ  34,5, END=3  Z5  I 
BEAD  35,5,END  =  3  Z6  I 
BEAD  36,5, END=3  Z7I 
BEAD  3715,  END  =  3  Z8  i; 

5  FCEMAT  (F20.6) 

C  C  CNVEESICN  FROM  IN-CZ/AMP  TO  HEBERS: 


Z2  (I) 
Z3  I 
Z4  I 
Z5  T 
Z6  I 
Z7  I 
GC  Tl 


Z2  (I) *2. 54*9. 8/ (100. *16. *2. 2) *16.8  2/5. 1393 
Z3  I  *2.54*9.8/  100. *16. *2. 2) *16.82/5. 1393 
Z4  TV *2.54*9.  8/ (100. *16.  *2. 2) *16.  82/5.  1393 
Z5  I  /6C. 

Z6  I)/60. 

Z7  I  /€0. 


CONTINUE 
CAII  TEK618 

CALL  SHEBPA  ( *  MACMILLA  '  ,  » A  *  ,3) 

CAII  PAGE  (8.5,11.0) 

CALL  NOBRDR 

CAII  PHYSOB  (1  .  ,  5.5) 

CAII  ABEA2D  (  5. 0,2.0) 

CAXL  COMP1X 

CAII  BEADIN  {'  MAGNETIC  FLUX  MODEL! » , 100 , 1. 5 , 2) 
CAII  HEADIN  {'ERASE  A  FLUX$ •  .  1 00,  1 .  ,  2) 

CAII  XNAME  PP BASE  B  FIOX$',jOO) 

CAII  YNAME  ('  $*,100) 

CAII  XIHTAX 

CAII  GBAF  (0.0, 30. ,360. ,-0.10, 'SCALE', 0.10) 
CAII  COHVE  (Z1,Z2, 629,0) 

CAII  DOT 


CAII  COBVE  (Z1.Z5, 629,0) 

CAU  BESET  (*  DOT') 

CAII  ENDGRjO) 

CAII  PHYSOB  p..  3.25) 

CAII  ABEA2D  (  5. 0,2.0) 

CAII  CCHPIX 

CAII  XNANE  (* PBASE  C  FLUX$*,100) 

CAII  YNAME  ( *  F1DX,  HEBERSl' , 100) 

CALL  XIHTAX 

CAII  GBAF  (0. 0,30.. 360. ,-0.  10, 'SCALE', 0.10) 
CAII  CURVE  (Z 1 , 23,529,0) 

CAII  DOT  _ _ „ 

CAII  COBVE  (Z1.Z6, 629,0) 

CAII  BESET  ('DOT ') 

CAII  ENDGRJ01 
CAII  PHYSOB  (1.  ,  1.) 

CAII  ABEA2D  (5. 0,2.0) 

CAII  XINTAX 

CAII  XNAME  ('BCTATION  ANGLE,  DEGREES!' ,  100) 
CAII  YNAME  J*  I'.IOOf  „ 

CAII  GBAF  (0. 0,30. .360. .-0. 10, 'SCALE' ,0. 10) 
CAII  CURVE  (Z1,Z4,629  ,0) 

CAII  DOT 

CALL  COBVE  (Z 1 .27 ,629, 0) 

CALL  RESETPDOI') 

CAII  ENDPL(O) 

CAII  CONEPL 

STCP 

END 


(f  *,100) 

;$* ,  loo) 
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THIS  PROCEDURE  SERVES  SEVERAL  PORPOSES.  FIRST,  IT  PEOVIDES  THE  MODEL 
MEANS  OF  APPLYING  THE  P HOPES  VOLTAGE  TO  THE  APPROPRIATE  WINDINGS  WHICH 
SIMULATES  COUNTERCLOCKWISE  COMMUTATION. 
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